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' For the Stackel family of the integrable systems a non-canonical transformation of the time 

variable is considered. This transformation may be associated to the ambiguity of the Abel map on 
the corresponding hyperelliptic curve. For some Stackel's systems with two degrees of freedom the 
2x2 Lax representations and the dynamical r-matrix algebras are constructed. As an examples 
the Henon-Heiles systems, integrable Holt potentials and the integrable deformations of the Kepler 
problem are discussed in detail. 
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In her celebrated papers |Q] published in 1889 S.Kowalewski discovered a necessary condition 
for an n-dimcnsional system to be completely integrable. This criterion enabled her to classify 
all integrable solid body motion about a fixed point and introduce the separation variables for 
her celebrated top. Evidently, it was the first application of singularity analysis to a concrete 
physical problem. Recall, that the method of singularity analysis associates integrability with the 
Kowalewski-Painleve property, i.e. a movable polelike singularity (t — io) _m in the solution of the 
equations of motion, here space and time must be thought of as complex [|. 

There exist cases, however, of integrable hamiltonian system with rational integrals of the 
motion, whose analytic structure permits solutions with algebraic singularities of the type (t — 
to) 1 ' 1 *, (k being a positive integer larger than one). This led to the introduction of the "weak" 
Painleve property A simple change of the independent variable, in and of itself, does not 
turn a "weak" Painleve system into one that satisfies the usual Kowalewski-Painleve criterion. In 
particular, the simple-minded idea to take (t ~ to) 1 ' 1 ' as a new independent variable does not lead 
to a Painleve expansion for all the solutions. According by any transformation that modifies 
the nature of the singular expansions must also involve a change of the dependent variables in order 
to reestablish the Painleve property. Such transformations, however, if they exist, are expected to 
be quite nontrivial and difficult to generalize to other examples. 

Of course, transformations, which use change of the independent time variable t, are non- 
canonical transformations. The classical example of the such non-canonical transformations is the 
duality of the Kepler problem to the geodesic motion on the sphere [Q or to the harmonic oscillator 

§■ 

Another important to us example is the Kolosoff transformations for the Kowalewski top 
@. In this case the Kowalewski separation variables coincide to the standard elliptic coordinates 
{qi > 92} at the plane (x, y) after the non-canonical change of the time 

dt= (q 1 (t)+q 2 (t)}dt. (1.1) 

Recently, the non-canonical transformations relate the Kowalewski top with the geodesic motion 
on SO (4) Q and with the Neumann system on the sphere S 2 

In this paper we consider the fortiori integrable systems with the known separation variables. 
Since we shall not discuss singularity analysis j^] and theory of algebraic completely integrable 
systems [pi in detail. Recall, the variable separation method permits one to reduce an integration 
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problem with several degrees of freedom to a sequence of one-dimensional integration problems. 
The inverse problem of obtaining various classes of completely integrable systems starting from 
a set of separated one-dimensional problems was started in the lectures by Jacobi. In framework 
of this approach, Liouville and Stackel introduced a family of the simple systems integrable in 
quadratures (a Liouville family is a particular case of a Stackel family). 

Here for the Stackel system we introduce non-canonical transformations of the time variable 
associated to the ambiguity of the Abel map on the hyperelliptic curve. All these transformations 
depend on coordinates only and, therefore, they are closed to the Kolosoff change of the time 
(1.1). For the some Stackel systems we propose the Lax pairs and r-matrix algebras. As an 
examples the Hcnon-Heiles systems, integrable Holt potentials and the integrable deformations of 
the Kepler problem are discussed in detail. 



2 Duality between the Stackel systems 

Before proceeding father it is useful to recall the classical work of Stackel || . The system associated 
with the name of Stackel Q is a holonomic system on the phase space R ™ equipped with the 
canonical variables {pj, with the standard symplectic structure f2„ and with the following 

Poisson brackets 

n 

Cl n = ^2 d Pj A d 1j j {Pj,Qk} = 5j k ■ (2.1) 

3=1 

The nondegenerate n x n Stackel matrix S, whose j column Skj depends only on qj 

detS^O, 7^ = 0, j^m 
oq m 

defines n functionally independent integrals of motion 

n q 

h = J2 c i« (P* + U i) > c ^ = dels ' (2 - 2) 

3=1 

which are quadratic in momenta. Here C = [cik] denotes inverse matrix to S and Skj be cofactor 
of the element . 

Each integral Ik ( [2.2| ) may be associated to the time variable tf., such that for any function 
£(p,q) one gets 

—77 = Ufc,£(p,q)}- 

dt k 

By definition the first integral I\ — H be the Hamilton function associated to the time t. The 
common level surface of the integrals fl2.2| ) 

M a = {z G M 2 " : h(z) = tt fc , k = 1,. ..,n} (2.3) 
is diffeomorphic to the n-dimensional real torus and one immediately gets 

where S{q\ . . . , q n ) is a reduced action function [ fTo] . If this real torus is a part of complex algebraic 
torus, then the corresponding mechanical system is called an algebraic completely integrable system 
§• 

The corresponding Hamilton- Jacobi equation on M a 



OS TTI dS OS . n dS dS 

-7W + H ^B7'---'B^' qi '---' qn > = ' ^ Cll 1^7ur 

oz oqi uqn ®Qj @Qj 

admits the variable separation 



S(qi ...,q n )=^2S j (q j ), Sj(qj) = / ^Fj(qj) dqj . (2.6) 
3=1 J 
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Here the functions Fj(X) depend on n parameters {afc}^ =1 

n 

F j (\) = Y / <* k s kj (\)-U j (\). 



k=l 



Coordinates qj(t, a±, . . . , a n ) are determined from the equation explicitly depending on time 



E 



s\ 3 (X)dX 



3 jt-/70(P0,3o) VJ2k=l a kSlj(X) ~ Uj{\) 

and from other n — 1 equations 



E 



^•>7o(po,«a) OfcSfei(A) - C/j-(A) 



(2.7) 



(2.8) 



The solutions of the problem is thus reduced to solving a sequence of one-dimensional problems, 
which is the essence of the method of separation of variables. 

Now we turn to the non-canonical change of the time and prove the following 

Proposition 1 // the two Stdckel matrices S and S be distinguished the first row only 

s kj = s kj i k 7^ 1 , 

the corresponding Stdckel systems with the following Hamilton functions 

detS(gi, . . . ,q n ) 



H = v(q)H, v(q) 



are related by non-canonical change of the time. 



det S(gi, 



(2.9) 



In fact, the corresponding Hamilton functions H and H obey to the equation (|2.9|), which follows 
from the definitions of the hamiltonians 



(2.10) 



and entries of the inverse matrix 



Cji 



1 9detS 



detS detS ds 



Equation (2.S) defines an implicit change of the time t —>■ t associated to the integrals H = 
I\ and H = Ii, respectively. On the other hand the equation (2.7) may be considered as an 
explicit determination of this transformation t — > t. In contrast with the general coupling constant 
metamorphosis discussed in J2| equation (2.9) is independent on the any constant entering in the 
potential U. 

Obviously, by using row by row transformations of the Stackel matrices with the associated 
tk — * tk transformations we can reduce the given Stackel system to any other Stackel system on 

HT)2n 



As an example, let us consider three matrices 



1 1 

1 -1 



<7i 92 
1 -1 



). i-l 


/ q\ 


i) 









(2.11) 



The corresponding hamiltonians H, H and H defined by ( 2.1CQ are dual (2.9) 



H 



H 



H. 



(gl+g|) 1 TT gl +g2 5 
£1 — t 7J tl . 

2 q\ + q\ 



(2.12) 
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for any potentials U. For t he h amiltonians H and H the change of the time ( 2.12 ) is closed to 



the Kolosoff transformation (1.1) |6[ and for any function £(q) depending on coordinates only one 
gets 

dt 2 2(qi+q 2 ) at 

For instance, let us consider uniform cubic potential 

Ufa) = 2a 2 q) +Pq] + ~f qj +6, (2.14) 

which gives rise to the hamiltonian H 

h = \(jp\ + vl) + a 2 (qf + ql) + ~ (q 2 i + <&) + \{<a + i2) + 5. (2.15) 

Using canonical transformation 

x + y 

91 = — — , Pi =Px+ Py , 

x ~~ y 

92 = — — , P2 =Px-Py, 

for the first system, the more complicated transformation 

31 = 7 x2/3 + P- > pi =Px x 1/3 - v ' 

4 d a 2 

92 = 7 z 2/3 -f^, p 2 = x 1/3 + y , 
4 o a 2 

for the system associated to S and the following change of variables in the third case 



Qi = V% - Vv , pi =PxVx-p v Vy, 

92 = -i (s/x + y/y) , p 2 — i (p x Vx + Py Vx) . 



one gets the Hamilton functions in the natural form 



H = ^ipl+pD + ^xix^+S^ + ^x' + y^ + lx + .S. 



H = 1 -( P l +P l) + 9 -^ x -y3( 3 - x 2 +y i) + S X -^ + ?l 1 onlyby ft = , 
# = o^Py-Tr— == + 7 — : 7T +<5 ' b y " = 0. 



2 ^ 4 V Vv 

under restriction (3 — for the second case. 

The system with the first Hamiltonian H is so-called first integrable case of the Henon-Heiles 
system [0. The second Hamiltonian H is related to so-called Holt potential J2j. Note, the second 
integral of motion is the polynomial of the third order in momenta for the Holt system. The 
system with the third Hamiltonian H may be considered as an integrable deformation of the 
Kepler problem. 

Duality between the Henon-Heiles and the Holt systems with the Hamiltonians H and H 
may be considered as the known coupling constant metamorphosis with respect to the constant 7 
||. The second known duality between the harmonic oscillator and the Kepler problem with the 
Hamiltonians H by a — ,7 = and H may be considered as the coupling constant metamorphosis 
with respect to another constant 8 [|[ |) . 

We can see that in practical circumstances the Stackel approach is not very useful because 
it is usually unknown which canonical transformation have to be used in order to transform a 
Hamiltonian ( 2.1C ) to the natural form H = T + V . This problem was partially solved for the 
uniform systems Uj — U, j = 1, . . . , n with polynomial potentials by using the corresponding Lax 
pairs JO]. Note, that the movable branch points of the type (t — t^) 1 ^ appear in the expansions 
of the physical variables [x, y) after canonical transformations. 

Henceforth, we shall restrict our attention to the uniform Stackel systems, where all the 
polynomial potentials Ujfa) = Ufa) and associated hyperelliptic curves Cj (2.4) are equal. 



4 



3 Duality and Abel map. 

Let us briefly recall some necessary facts about the Abel map and the inverse Jacobi problem. 
The set of point C(z, A) satisfying 

2g+l 2g+l 

C: z 2 — F(X) — GuX k — (A — Xj) , (3.1) 

k=0 j=l 

is a model of a plane hyperelliptic curve of genus g. Here F(X) is polynomial without multiple 
zeros. Let us denote by Div(C) the Abelian divisor group and denote by J(C) the Jacobian of the 
curve C. The Abel map puts into correspondence the point D € Div(C) and the point u £ J(C) 

U : Div(C) -> J(C) , (3.2) 

The set of all effective divisors D — 71 + • • • + 7 n (the may be not mutually distinct) of degn 
of C is called the nth symmetric product of C, and is denoted by = S n C. The C^ n > can be 
identified with the set of all unordered n-tuplcs {71, . . . ,7 n }, where jj are arbitrary elements of C. 



Now consider restriction of the Abel map (3.5) to 

U : C {n) -» J(C) , (3.3) 

where 

W(7i,7a , • • • ,7 9 ) = W(ti) + U{ l2 ) + ■ ■ ■ +U( lg ) . 

According to the Abel- Jacobi theorem this map is surjective and generically injective if n — g only 
[jjj], |l3). If n g the Abel map is either lack of uniqueness or degenerate. The corresponding 
Stackel system either has a dual system associated with the same curve or it is a superintegrable 
system [p"l|| . 

S upp ose that point D = 71 + • • • + jk, k < g belongs to C^ k ' . The differential of the Abel- Jacobi 
map ( |3.3| ) at the point D is a linear mapping from the tangent space Td(C^) of at the point 
D into the tangent space T U ^(J(C)) of J(C) at the point U(D) 

W D : T D {C^)^T U(D) {J{C)). 

Now suppose that D is a generic divisor, and Xj is a local coordinate on C near the point jj. Then 
(ari, . . . , x n ) yields a local coordinate system near the point D in C^ n \ Let dw^ (k — 1, ■ ■■,<?) is a 
basis for a space Ti.i(C) of holomorphic differentials on C, and near 7^ 

dw fe = 4> k3 (xj)dxj , (3.4) 

where <fikj{xj) is holomorphic. It follows that the Abel- Jacobi map Li can be expressed near D as 

Hence 

(0n (71) ••• ^1(71) \ 

: •.. : . (3.5) 

<t>lk(7n) ■ ■ ■ (Pgniln) J 

is the so-called Brill- Noether matrix Jl4| . Hencefort h, w e shall restrict our attention to the special 



divisors D s , such that coefficients in the expansion (3.4) are independent on the point jj 



dw k = (f> k (xj)dxj . 

In this case all rows of the symmetric Brill- Noether matrix depend on local coordinate {x\ , . . . , x n } 
identically. 



The Jacobi inversion problem (2.8)) is formulated as follows: for a given point 



u=(f3 1 ,f3 2 ,...,f3 n ) e J(C) 
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find n points 71 , 72 



, 7„ on the genus g Ricmann surface C such that 
a 



E 



dw-j 



j = 1, .. . ,n. 



(3.6) 



fe=l 17 T° 



Here we shall tacitly assume that the base point 70 € C has already been fixed 12 . 

If n = g for almost all points u S J(C) the solution D = 71 H + 7„ exist and is uniquely 



determined by system fl3.q ) (for the unordered set of points jj) ]12|]. However, if the degree n < g 
of the symmetric product C^ n ' is less than genus g of C, the Abel map is lack of uniqueness. 
In this case we can propose that two different points u, u € J(C) have the one Abel preimage 
{71,..., 7n} 6CW. 

The Abel preimage of the point u e J(C) is given by set {{pi, gi), . . . , (p n , q n )} & C^ n \ where 
{qi, . . . , q n } are zeros of the Bolza equation |ll], [l| 



e(A,u) = A n - A"- 1 P n,n(u) - A n - 2 P «,n-l(u) 

and {pi, . . . ,p n } are equal to 

Pk - 



p„,i(u) = 0, 



de(\, u) 



A=«Jfc 



Here vector u belongs to Jacobian J(C) and pkj(u) is the Kleinian p-function [|T5|, |H 
Now we turn to the uniform Stackel systems. We can regard each expression (p 
defined on the genus g Riemann surface 



C 



F(X) 



(3.7) 
(3.8) 
as being 

(3.9) 



fc=i 



which depends on the values <Xk of integrals of motion. For the Stackel systems on R the 
minimum admissible genus g of the curve C is equal to g = [(n — l)/2]. 

The nth symmetric product of C defines the n-dimensional Lagrangian submanifold in the 
complete symplectic manifold M. n 



C (n) 



C{pi,qi) X C(p 2 , 92) X • • • X C(j>n, qn) ■ 



(3.10) 



Then, the integration problem (2.7 -pT^ ) for equation of motion is reduced to inverse Jacobi problem 
(3.3) on Lagrangian submanifold ( |3.10 ). The corresponding holomorphic differentials <5w/fc are equal 
to 

SkMdX (3U) 



dwk 



z(A) 



The set of these differentials either form a basis in the space of holomorphic differentials 7i\{C) 
Ji"2| or may be complement to a basis. The corresponding n x n Stackel matrix be the n x n block 
of the transpose Brill-Noether matrix . 

The different blocks are determined the dual Stackel systems. In this case vectors differing 
the first entry only 

u = {t, fo, . . . , p n } £ J(C), vi=fclh,... > p n } e J(C) 

have a common Abel preimage {(pi, • • • , (Pn,q n )} € C^ 71 '. 

Let us consider the standard basis of holomorphic differentials in Hi(C) 



AJ- 1 
aw,- = —7-t- dX . 
3 z(\) 



J = 1, 



,9- 



(3.12) 



Recall, that derivative Up bears a great resemblance to the canonical map C — > P 9 1 and, therefore, 
to the Veronese map given by a basis for the polynomial ring of degree g — 1. With 

respect to the basis of H%(C) ( 3.12 ), the Veronese map of C has an extremely simple expression 

(y,\)^ \^l\°-\\*- 2 ...,A,1]. 
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By using the corresponding symmetric Brill-Noether matrix Up (f3Jj), we shall determine the 
Stackel matrices as (n x n) blocks of the following (g x n) matrix 



/ 9? 



9-1 



9-1 
<?2 



n-2 (7-2 
<7l <f 2 



V i 



-,3-1 



7 9-2 



(3.13) 



1 / 



Evidently, all the Stackel matrices can not be obtained from the symmetric Brill-Noether matrices. 
For instance, the Stackel matrices (2.11) do not belongs to the set of symmetric matrices. 



4 Lax representation. 



Henceforth, we shall restrict our attention to the basis ( 3.12 ) and the symmetric matrix ( 3.13 ). 
For the corresponding Stackel systems let us look for the Lax representation as 



L = 



/i(A,p,q) e(A,q) 
/(A,p,q) ~h(X,p,q) 



(4.1) 



Hereafter, by abuse of notation, we shall omit the some arguments at the entries of the Lax matrix. 

Let us fix hyperelliptic genus g curve C and dimension of the phase space n < g. Then we 
extract the (n x n) Stackel matrix S from the matrix ( 3.13| ) and define the Hamilton function H 
( |2.10D with [7 = 0. 

To construct the Lax matrix let us determine function e(A, u) ([3.7]) initially 



e (A,q) = ;Q( A -^), 
3=1 

with n zeroes, which are solution of the inverse Jacobi problem. 

In the second step let us introduce the second entry of the Lax matrix as 



(4.2) 



h(X) 



1 



de(X) 



+ w(A,p,q)e(A). 



2«(A,q) dt 

Here function v(X, q) is calculated by using the second Bolza equation ([sl| 



h(X)\ x =Pfc = (l^h 
v ;i x= qk Fk y 2v dt j 



de(X) 



2v dt 

Let the third entry of the Lax matrix takes the form 



A=9fc 



(4.3) 



A=9fc 



v dt 

Here the single unknown function w(X, p, q) is determined such, that the spectral curve of the Lax 
matrix (4.1) 

C: z 2 — F(X) — — det L (X) — h 2 (X) + e(X) /(A) (4.4) 



be the same as initial algebraic curve C (2.4) by U = 0. 
The constructed above matrix Lg(X) (4.1) reads as 



L (X) = 



( -i-e t (\) +w(A,p,q)e(A) 
2v 



where 



de(X) 
dt 



h t (X) 



= {H,e(X)}. 



1 



h,= 



e(A) \ 
e t (A) - w(A,p,q)e(A) / 



(4.5) 
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obeys the Lax equation 



with the second matrix 



Aq, Lq 



v(X,q) 



™(A,p,q) 1 

-w(A,p,q) 



By definition of the Lax matrix all the pair s of separation variables jj — (pj ,qj) (4.2- Ljj ) lie on 
the spectral curve C (4.4) of the matrix (|4.5|) 



z 2 ( 7j ) = p= = ^ 2 (A)| A= ^ = F(A = 5i ) = F(A)| 7j . 



For the systems with polynomial potential U ^ we propose to change the entry /(A) in (L5) as 

/(A) = i^+u(A,q)e(A), 

where we add new function u(A, q) depending on coordinates only. Of course, to construct the Lax 
matrix here 



L(X) = 



( -7T e*(A) + io(A, p, q) e(A) 



e(A) 



(4.6) 



V - /i t (A) +u(A,q)e(A) — e t (A) + iu(A,p,q) e(A) / 

we have to use the complete Hamiltonian with U ^ 0. The associated second Lax matrix reads as 



A= An 





v(X,q)u(X,q) 



w(A,q) 



w(A,p,q) 1 
u(A,q) -u)(A,p,q) 



(4.7) 



To consider the corresponding Lax equation 

dL(X) 
dt ~ 



A(X),L(\) 



we can assume that the common factor w(A, q) in front of the matrix A may be associated to the 
change of the time for the Stackel systems. 

In general the proof of existence functions v, w and u require s an applicatio n of the method of 
algebraic geometry [13]. By definition of the Lax matrices L(X) (4.6) and A(X) (LI ) this problem 
may be reduced to the solution of the single equation 



df(X) 
dt 



2v(uh-wf) =0, 



dF(X, e, v, u) 
dt 



0, 



(4.8) 



for the given function e(A) ( |4.2| ) and the given Hamiltonian H ( |2.10 ). 

If we consider the lower (n x n) block of the matrix (3.13), the differentials ( 3.1 1| ) span a whole 
space TL\(C) and the Abel map is the one-to-one correspondence. In this case from equations (3.8) 
and (4.3) follows that 

v t (X,q) = 0, w(A,p,q)=0. 



If we put v = 1, rename t — x and introduce "new" time variable r, the equation (4.8) is rewritten 



du(x, t, A) 



ld 3 x + u(X)d x + ~u x (X) 



e(A) = , x = t . 



(4.9) 



This equation may be identified with equation on the finite-band stationary solutions = 

of the nonlinear soliton equations. In this theory equation (4.9) is called the generating equation. 
For different choices of the form of e(A) and u(X), this procedure leads to different hierarchies of 
integrable equations, as an example to the KdV, nonlinear Shrodinger and sine-Gordon hierarchies 
or to the Dym hierarchy (see references within JTT[| ) . 

Function u(X, q) in ( |4.9| ) is constructed by using function e(A) (|3jJ|4]| 



u(X,qi, . . . ,q n ) = (j)(X)e 2 (A) 



MN 



(4.10) 
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Here 0(A) is a parametric function on spectral parameter and [£]jv is the linear combinations of 
the following Taylor projections 



+00 



JV 



N 



J2 z ^ k 



(4.11) 



N 



k=0 



or the Laurent projectio ns |p7| , [Tl| . 

If the differentials ( 3.11 ) span the whole space Hi(C) the corresponding Stackel systems de- 
scribe all the possible systems, which separable in the orthogonal curvilinear coordinate systems in 
M." pHl . Let us consider the Stackel systems which are d ual to these systems. To apply equation 
( 2.13 ) to the function e(A) ( fl.2| ) and by using definition (L3) one gets 



p k = h(X) 



-2~{^e(A)} 



detS 
detS 



(-> 



Pk 



detS 
detS 



detS 
detS 



(4.12) 



Recall that v = 1 for the integrable system with the Hamiltonian H associated to the lower {n x n) 

block of the matrix ( 3.13 ). 

Thus, according to ( 4.12 ), below we shall consider the Stackel systems with following functions 
v(q) only 

detS(gi, . . . ,q n ) 

v{q) = . 

det S(gi, . . . , q n ) 



The correspondi ng c hange of the time ( 2.13 ) depending on coordinates only is closed to the Kolosoff 
transformation (IT) |J. 

Let us briefly discuss canonical transformation which transforms a Hamiltonian (2.10) to the 
natural form H = T +V . For integrable systems separable in the orthogonal curvilinear coordinate 
systems on M. the Abel map is one-to-one correspondence and vt = {H, v} = 0. In this case we 
can put v = 1 and introduce function 0(A) 



# 2 (A)=e(A), 



(4.13) 



which was proposed in the theory of the soliton equations [[16| . It allows us to rewrite generating 
function of integrals of motion 



F(X) = — £> 3 Bu + u(X, q) B A . 
as a Newton equation for the function B 

B(X, q) = —F(X, a u ...,a n ) B- 3 (X, q) + u{X, q) B(X, q) 
To expand function B(X) at the Laurent set 

N 



(4.14) 



(4.15) 



0=0 



XjX J 



it is easy to prove that coefficients Xj obey the Newton equation of motion (4.15) (see references 
within |l(], |ll|). Here we reinterpret the coefficients of the function F(X) in (4.15) not as functions 
on the phase space, but rather as integration constants ctj (2.3). 

In general by v t ^ the generating function F(X) = — det L(X) (4.6) is equal to 



F(\) 



4 M^- 2 "''M^-«')v+(^) 



In this case the suitable canonical transformations, which transforms any Hamiltonian (2.10) to 
the natural form, are unknown. 

Although we can not proof validity of the presented Lax representation in general, this con- 
struction works for the many well-known mechanical systems. In the next Section we consider 
some two-dimensional Stackel systems in detail. 
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5 Examples. 



Let us consider four orthogonal systems of coordinates on plane: elliptic, parabolic, polar and 
cartesian |^. The Lax matrix Lq(\) (4.5) by U — is transformed to the Lax matrix L(X) (l.t) 
by U ^ by using the outer automorphism of the space of infinite-dimensional representations of 
underlying algebra sl(2) jl?], [H]. Since, we shall consider the Lax representations for the geodesic 
motion by U — more extensively. 

1. Parabolic and cartesian coordinate systems (w(X,p, q) = 0). 
Let us consider two hyperelliptic curves 



C( 2 )2 



2g+l 
i=l 

2g+l 

=x- 1 H(\-\ t ). 



(5.1) 



If we choose the standard basis in the space of holomorphic differentials one gets the following 
symmetric matrices ( 3.1 3| ) for two-dimensional systems 



U\ (91,92) = 



9? 



4^ \ 



4 



9i 92 

V -1 -1 / 



(91,92) = 



/ 9-2 

/ 9i 



9i 
1 

1 

9i 



V 



9f- 2 \ 



92 

1 

1 

92 ) 



(5.2) 



Different (2 x 2) blocks of the matrices U* 1 determine different Stackel systems. 

Let us consider two blocks for the each matrices, such that the corresponding change of the 
time will be same as the Kolosoff transformation (1.1) Q. So, for the curve we shall consider 
the following matrices 



Si = 



9i 92 
-1 -1 



Si = 



? 9 

9i 92 



-1 -1 

For the second curve the associated Stackel matrices are equal to 

! 1 \ / 91 92 



(5.3) 



S. - , ^ ^ 

9i 92 

Introduce the Hamilton functions ( 2.10[ ) by U = 



V 



1 

91 



1 

92 



(5.4) 



H, 



(i) 



H, 



(2) 



2 2 

H -P2 
9i - 92 



9i Pi - 92 v\ 



^ 1J =(9i+92)- 1 ^ 1 



(5.5) 



(2) 



9i - 92 

The corresponding second integrals of motion of the dual systems are related 

m _ T (k) 



J Q — Jq 



9i 92 „(fc) 

— ; — -"0 ' 

9i+92 



The functions e(A, u) (|J) 

ei(A) = (A-9i)(A-9 2 ), 



e 2 (A) 



A: = 1,2. 



(A - giKA - q 2 ) 
A 



(5.6) 
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have two zeroes, which are solution of the inverse Jacobi problem ( pTg| ) on C*- 1 -* and , respe ctive ly. 
Let us introduce the new physical variables at once. For the first curve C' 1 ' equation ( 4.15| ) 



ei(A) = (\-qi)(\-q 2 ) = B 2 (\) , B{\) = A 
immediately yields the following canonical transformation 



^__y_ 

2 4A 



q-2 



x + y/2y 



Pi = Px 
P2 = Px 



2yp y , 
Zypy ■ 



These variables obviously related to the cartesian coordinate system. For the second curve 
the corresponding equation 

e 2 (A) = A" 1 (A - gi )(A - <? 2 ) = A - x - ?L 



defines the standard parabolic coordinate system 

x - 



'72 



\/x 2 + y 2 



Pi = Px 



x 2 + y 2 + . 

y 



■Pv. 



x + \J x 2 + y 2 



By [7 = the Hamilton functions are given by 



H, 



(i) 



\/ x 2 + y 2 — x 

P2=Px-\ Py ■ 

y 



(2) 



'0 ~^PxPy, H(f' =Px+Py- 

According to (4.3) and ( [1.12| ) functions v{q%, (72) entering in the Lax representation are equal to 
u (A, (ft, (72) == 1 for matrices Si t 2 , 



(5.7) 



v[\qi,q2) = (qi +92) 1 = 
In physical variables the Lax matrices are given by 



1 



for matrices Si. 2 • 



-(i) 



p x + (2A — x)p y A 2 — A x + 



x 2 - 2 y 



-p x - (2A - x) p y 



/ 1 x 1 2 \ 

/ Px + T^yPy A - x- —y z \ 



(5.8) 



-(2) 



2A 



4A 



V 



1 



1 



~Px 



x ry rx Py J 

— (1 2) 

For the dual Stackel systems the Lax matrices L y ' ; have the form 





L 



(i) 



.PxPy 



1 




( ° 










:) 


0) 





(5.9) 



/ \ 



L, 



(2) _ t(2) 



V 



PS + 















By using property {h t (X),v(q}} = of the function v(q) (|5.7[) we can easy proof equation ( |4.8| ) for 
the dual systems by using the same equation for the system with v t = 



f t = {H,f + H} = ft = 0. 



11 



Another consequence of this property is that the function w(A,p,q) in (4.5) is equal to zero. 

Note, in the works and S, devoted to the Kowalewski top, the common Lax matrices are 
proposed for the both dual systems after the non-canonical change of variables. Here we obtain 
different Lax matrices for the systems connected by non-canonical change of the time. 

i 1 f-t 2) j 1 

The spectral curves of the matrices Lq (|5.8| ) coincides with the initial curves Cq ' (|2.4| ) by 
U = 



(i) 



z A = H, 



(2) 



Ji 



(2) 


A 



(5.10) 



Here Jg 1 ' 2 "* be the second integrals of motion ( |2.2| ). For the dual systems with the Hamilton 

~ (1 2) 

functions ' the corresponding spectral curves are equal to 



z 2 = A 2 + J, 



(i) 
o > 



z 2 = H^X 



7(2) 

'o 
A 



(5.11) 



If for the system with the Hamiltonian Hq the Abel map is one-to-one correspondence on the 
curve ( 5.10| ), then for the same system on the curve 

z 2 = e 2 A 2 + ei A + e 

the associated Abel map is lack of uniqueness in general. So, on this curve we can introduce the 
second Stackel system with the dual Hamiltonian Hq~\ 

Let us briefly consider systems with polynomial potentials U 7^). As an example, introduce 
different potentials for the curves C^ 1 ' 2 ) (5.1) 



9j 



(5.12) 



To describe these potentials we have to put N = 6 and N = 4 in ( 4.11 ) and have to use the 
following parametric functions 

0(D (A) = -a 2 A 5 and <^ 2 )(A) = -a 2 A 3 

for the curves and C^ 2 \ respectively. For the both curves the common function u{\, 51,(72) is 
given by 

u {1 ' 2) =~a 2 {\ + 2x). (5.13) 

Here we restrict ourselves the pr esen tation of the function u only, the complete Lax matrices L(X) 
may be constructed by the rule (4.6). 

The spectral curves of the corresponding matrices ([0]) coincides with the initial curves ( |3.9| ) . 
For instance, curves for the systems with dual Hamiltonians H 1 -- 1 ' 2 ^ are 



C« : 
: 



z = a \ J + {3 \ — H \ — J . 



a 2 A 3 -HX + /3 



J 



The Poisson bracket relations for the Lax matrix (|5~8-5.9) are closed into the following linear 
r-matrix algebra 



{L(A),L( M )} = [r 12 (X,fi),L(X)} - [r 21 (A, M ),2( M )], 



(5.14) 



^2i(A, (i) = -nr 2 i(A,^)n. 

Here the standard notations are introduced: 



L(X) = L(X) <E> I , L(n) = I <g> L{n) , 

and n is the permutation operator of auxiliary spaces []l8[ . 

By Vt — for the systems related to the matrices S ' 2 the corresponding r-matrices ry(A,/Lt) 
in (5.14) consist of two terms 

nj=r^+rV.. (5.15) 
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The first matrix is a standard r-matrix on the loop algebra C(sl(2)) 



n 



A — fi A — ji 



( 1 \ 

10 

10 

V o o o i l 



(5.16) 



The second matrix may be associated to outer automorphism of the space of infinite-dimensional 
representations of underlying algebra sl(2) |T^, [H]. The corresponding dynamical r^-matrices 
depend on the coordinates only 



u(X,q) - u(fj,,q) 
r i2 = : 





1 



(5.17) 



A — /! 

By v t 7^ for the dual Stackel systems related to the matrices S12 we have to add to the 



r-matrices 5.15 the third term 



' 12 



( 











\ 

















1 










V 











) 



(5.18) 



where the second matrix r\ x is defined by ( |5.14 ) . 

The matrix t\a may be connected with the Drinfeld twist for the Toda lattice associated to 
the root system T> n . Let us consider the Drinfeld twist ]19[ of the quantum R- matrix 



R = FRF 21 1 , F 2 i=UFn. 



(5.19) 



Here matrix R satisfies the Yang-Baxter equation and matrix F has the special property J19| . To 
introduce the corresponding linear r-matrix poj , one gets 

R = 1 + 2r)r p + 0(rf) , F = I + rj r° + 0{r/ 2 ) . 

Then we consider limit of the twisted matrix R by r\ — > 

R 12 =I + v(r P 12 + r\ 2 - B ( r p 12 + r\ 2 ) II ) + 0(rj 2 ) . 

Formally, coefficients by 77 may be called twisted linear r-matrix. 
By using generators h, e, f of the underlying Lie algebra sl(2) 



(5.20) 



(5.21) 



[h,e]=2e, [h,f] = -2f, [e,f]=h, 
let us introduce an appropriate element T £ U(sl(2)) ® U(sl(2)) 

Ft = exp(t • e f) , C e C 

belonging to a tensor product of the corresponding universal enveloping algebras U(sl(2)) [jl9(| . In 
the fundamental spin-1/2 representation pi we have 



F(0 = (Pl®Pl)Fi 



To substitute in ( |5.19| ) the Yang solution of the Yang-Baxter equation R = 1+ — H we get a twisted 

A 

i?-matrix. If the element £(q) be a suitable function on coordinates, this dynamical twisted R- 
matrix may be used to description of the Toda l attice associated with the V n root system Q . 

Let us consider twisted dynamical matrix (5.19) by £ = f(q). We can see that the linear 
r-matrix associated to the dual Stackel system ( 5.16| - 5T§| ) 



/ 1 











\ 





1 













1 







V 








1 





ri2 



— ' 12 



12 



rai =-n(rf 2 



' 12 



)n 
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is equal to the half of the twisted linear matr ix ( 5.20| ). 

Recall, for the Stackel matrices Si (5.3) and S2 ( |5.4| ) the corresponding differentials (3.11) 
span Tii. Since, the associated Hamilton functions have a natural form in physical variables. For 
instance, Hamiltonians with potentials (5.12) are given by 

- 2 5 



%Px Py 



a 

T 



{y 2 + 5x 2 y+-x^ 



2 V 2 



■v) 



(5.22) 



P. 



V)+)9. 

To consider the dual Stackel systems we have to use additional transformation 

x \/2 x , Px — Px y/~2~x , 
for the first curve and the following more complicated transformation 



(5.23) 



y 



2py 

3 a 



Px Px % 



Py 



1/3 



ay 



for the second curve. After this canonical change of variables the Hamiltonians H^ 1,2 ) (5.5) obtain 
the natural form 



ftW = 2p xP y 



a 2 {y 2 + lOyx + 5x 2 ) 



P(3x + y) 



(5.24) 



{pi 



pI) 



3 QL — _0 /q ,9 — 

■ x A ' A (- x 



y 2 ) 



■/3x 



-2/3 



4 v 2 

The system with the Hamiltonian is so-called second integrable case of the Henon-Heiles 

system The dual system with the Hamiltonian is so-called Holt-type system Note, 
the second integral of motion is a polyno mial o f the fours order in momenta for the Holt system. 

Additional canonical transfor mati on ( 5.24 ) allows us to get natural Hamil tonia ns for the re- 
stricted class of the potentials U ( 5.12 ) only. Unlike canonical transformation ( 5.23] ) may be used 
for any potentials U. As an example, rational potential 



U(q): 

give rise the following Hamiltonian 

4a 



a 

7/2 



+ - +jq + Sq 2 + pq 4 



H = 2p x Py 



P 



(x - y) 2 



V 



y)- 



Also we can add potential terms (5.24) to this Hamiltonian. 

By v = 1 and w = the Lax representation ( |4.5| ) for a system with an arbitrary number n of 
degrees of freedom may be regarded as a generic point at the loop algebra C(sl(2)) in fundamental 
representation after an appropriate completion p~T| ] . As an example, for the generalized parabolic 
coordinate systems function e(A) is given by 



e(A) 



II"=i(A-g 3 -) 

n^(A-4) 



A 



n-l 



^4(A-4) 



6 k € 



To construct the Lax representation for a potential motion we can use the outer automorphism of 
the space of infinite-dimensional representations of sl(2) proposed in Jl7t . 

By Vt 7^ for the dual Stackel systems the Lax representations may be constructed without 
any problem as well. For instance, let us consider system with the three degrees of freedom. To 
construct the Lax matrix by (4.5-4.6) with the function u given by ( 5.1 3| ) one gets 



9 2 



H 



p 2 x+pI+pI + 



a k z 2 



k e . 

a 2 

+ Y 



2x 2 + y 2 + z 2 
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After an additional canonical transformation (5.24) extended on the p z , z variables the Hamilton 
function takes the form 



h = pi + f y + f z (i 



— 2/3 



3a 



r-2/3 



(o* +y + z ) 



So, the main unsolved problem is to introduce additional canonical transformation, which transform 
the dual Hamilton function H into the natural form. 

2. Elliptic and polar coordinates (w(p, q) 7^ 0). 

Recall, that the polar coordinate system may be obtained from elliptic coordinate system and, 
therefore, we shall consider elliptic coordinate systems in detail. 
For the elliptic coordinate systems algebraic curve is given by 



C (3) 



n^(A-A 4 ) 



(\-k)(\ + k) ' 
Let us consider two Stackel matrices associated to this curve 



fceC. 



/ 



'12 



s, = 



q\ — k 2 q\ — k 2 
1 1 



\ 



Sa = 



/ Aq 2 4q 2 \ 

q\ — k 2 q\ — k 2 



\q\- k 2 q 2 -k 2 J 



\ q 2 — k 2 q\ — k 2 J 



(5.25) 



The corresponding non-canonical change of the time (2.S) is closed to the Kolosoff transformation 
01 

For the polar coordinate system the Staackel matrices are non-symmetric matrices 



S 4 = 



/ if 



q\ A{q 2 -k) 



\ 



(5.26) 



\ q\ Hq 2 2 -k) J 



The corresponding non-canonical change of the time (2.£) is closed to the Kepler change of the 
time 

By U — the initial hyperelliptic curves (3.£) for the matrices S3 and S3 are given by 



H X + J 
\~ k 2 



z 2 = 



4ff A 2 + J 
A-fe 2 



(5.27) 



with the following Hamiltonians 



„(3) p{{q{-k 2 )~pi(qi-k 2 ) 
H Q = 



^ 3) = 



qi - qi 

The Hamiltonians related to the matrices S4 and S4 read as 

q\-k 



4(gi + g 2 ) 



H. 



Hi i] =pl-A. 



2 v\ 



H, 



(4) 



"2 ttW 



QiH { 
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(5.28) 



(5.29) 



Let us fix elliptic coordinates by using equation 



m _ (A - gi) (A - q 2 ) _ _ 
{ ' (A - k) (A + k) 



4(A-fc) 4(A + fc) 



such that 



Qi = 

92 = 



x 2 + y 2 



+ - y/(x 2 + y 2 ) 2 + 16 k (x 2 - y 2 ) + 64 k 2 



x 2 + y 2 1 



- ^(x 2 + y 2 ) 2 + 16 k (x 2 - y 2 ) + 64 k 2 . 
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The corresponding equation for the polar coordinates 



e(A) 



gi (A - g 2 ) 
A(A-l) 



x 2 Ay 2 
4A + A- 1 



immediately yields 



5i 



r = \Jx 2 + y 2 , <j 2 = cos 2 ((/>) 



In physical variables the Hamiltonians (5.28 - 572S| ) have a common form 



H - Px+Pyi H - 9 

x + y 



To cons truct the Lax representations we begin with the calculations of the functions v(X, q) 
by the rule (fD| fy]| 



for matrices S34 
1 



4 ' x 2 + 



(ft +92) 1 = 



for matrix S3. 4 



(5.30) 



q\ x 2 + y 2 



for matrix S4 



So, for the Stackel systems associated with the matrices S3 ( 5.250 and S 4 ( [5.26| ) one gets 



L (X) 



xp x 



VPy 



2 (A - k) 2 (A + Jfe) 

„2 



eA 



y 



\ 



P7, 



4 (A - k) 4 (A + jfe) 

XPx VPy 



\ 



(5.31) 



X-k X + k 2(X-k) 2(X + k) ) 

Here e = 1 for the elliptic coordinate system and e = for the parabolic coordinate system. The 
spectral curve of the Lax matrix Lq(X) coincides to the initial curve (5.27). 

For the dual system, in contrast to the cartesian and parabolic coordinates, the Lax matrix 
has the more complicated form. Both these Lax matrices may be constructed by the rule ( [4.6j ) 
with the following common function w(p, q) 



The Lax matrix reads as 



i (A) =L (A) + 



-2w 



w(p,q) = 2 VH 

w e(A) 
h(X) — w e(A) - eio 




e = 0,l. 



(5.32) 



(5.33) 



Here e(A) and h(X) are entries of the corresponding matrices L (X) ( p.31 ) by e = 0, 1. As above, 
the spectral curve of the Lax matrix L (X) by e = 1 coincides with the initial curve (5.27). 
For the cartesian and parabolic coordinate systems we can get equation 



{ { H, v-\q) } , e(A, q) } = { { H, (qi + q 2 ) } , e(A, q) } 



= 2. 



on the Hamiltonian H, function e(A) and function v(q) defining change of the time. For the polar 
and elliptic coordinate systems the corresponding equation is 



{{^, U - 1 (q)},e(A,q)}=8(e(A)-e) 



0,1. 



Hence, from the equation (4.8) follows that the function u>(p,q) in (4.5- 5.32 ) does not equal to 
zero. If we consider more complicated change of the time for the cartesian and parabolic coordinate 
systems, one gets non-zero function w (^^) as well. 



1G 



The Poisson bracke t rela tions for the Lax matri x £q (A) ( 5.31 ) are closed into the standard 
linear r-matrix algebra (5.14) with rational r-matrix (5.16|) on the loop algebra C(sl(2)) fll8f. 



The Poisson brackets relations for the dual Lax matrix Lq(X) ( |5.35| ) have a poly-linear form 



{L(X),L(p)} 



n 2 ,i(A) + r 21l L(y) 



+ RL(X) L(n) +L(X) L(ji) R -L(X) RL{^) -L{n)RL{X) 



Here linear r-matrix reads as 

r 12 (X^)=r p 12 (X-ri+4er? 2 , 



r2i(A,/i) = -nr 2 i(A,/j)n, 



where r p (X— fx) be the standard linear r-matrix on the loop algebra £(sZ(2)). The second dynamical 
term is given by 



( 











\ 














w 


1 










V o 


— w 











The quadratic R- matrix is closed to the twisted linear r-matrix (5.20) 



R = — [rf 2 
w 



2 

— i /•; 

w 



nrf 2 n 



For the systems with U ^ functions w(A, q) may be constructed as usual (jl7| Note, the 
both dual Hamiltonians obtain a natural form after the following additional canonical transforma- 
tion of variables 



x = vx-^/y, p x = Vxp x - \/yVy , 

y = -i (Vz + y/y) , p y =i {VSp x + y^p y ) 

As an example, for elliptic coordinate system the uniform potential 



give rise to the following dual Hamiltonian 



a 



H = 2 PxPy + -(x + k) (y + k) 







>xy 



(2 xy + k x + k y) . 



For the polar coordinate system we present the non-uniform degenerate potentials 



^ ( %i)=,3, 
associated to the dual Hamiltonians in the form 



Ut\q 2 ) = 0. 



H 



Px Py 



0( 



10 



1 1 

x y 



xy 



Both these systems may be considered as an integrable deformation of the Kepler problem. 



6 Conclusions 

In this paper we have considered the non-canonical relations between the different Stackel systems. 
The proposed change of the time is related to ambiguity of the Abel map. For the two degree of 
freedom systems that were studied in this paper, we found the Lax representations and the r- 
matrix algebras. The corresponding dynamical r-matrices have the intriguing connections to the 
Drinfeld twists. 



Of course, considered above particular family of the time transformations (2.9) does not ex- 
hausted all the set of the non-canonical changes of the time, which preserve the integrability. As an 
example, the complete Kolosoff transformation {t,p, q} — > {t,p, q} Q connects the Stackel system 
with the other integrable non-Stackel system. So, it would be interesting to investigate another 
integrable systems connected with the Stackel systems by non-canonical transformations. 

This work was partially supported by RFBR grant. 
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